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1.  SUMMARY 


A  seismic  moment  tensor  is  a  3  x  3  symmetric  matrix  that  provides  a  compact 
representation  of  a  seismic  source.  We  develop  an  algorithm  to  estimate  moment  tensors 
and  their  uncertainties  from  observed  seismic  data.  For  a  given  event,  the  algorithm 
performs  a  grid  search  over  the  six-dimensional  space  of  moment  tensors  by  generating 
synthetic  waveforms  for  each  moment  tensor  and  then  evaluating  a  misfit  function 
between  the  observed  and  synthetic  waveforms.  ‘The’  moment  tensor  Mo  for  the  event  is 
then  the  moment  tensor  with  minimum  misfit.  To  describe  the  uncertainty  associated  with 
Mo,  we  first  convert  the  misfit  function  to  a  probability  function.  The  uncertainty,  or 
rather  the  confidence,  is  then  given  by  the  ‘confidence  curve’  P(V),  where  P(V)  is  the 
probability  that  the  true  moment  tensor  for  the  event  lies  within  the  neighborhood  of  M 
that  has  fractional  volume  V.  The  area  under  the  confidence  curve  provides  a  single, 
abbreviated  ‘confidence  parameter’  for  Mo.  We  apply  the  method  to  data  from  events  in 
different  regions  and  tectonic  settings:  63  small  (Mw  <  2.5)  events  at  Uturuncu  volcano 
in  Bolivia,  21  moderate  (Mw  >  4)  earthquakes  in  the  southern  Alaska  subduction  zone, 
and  12  earthquakes  and  17  nuclear  explosions  at  the  Nevada  Test  Site.  Characterization 
of  moment  tensor  uncertainties  puts  us  in  better  position  to  discriminate  among  moment 
tensor  source  types  and  to  assign  physical  processes  to  the  events. 


2.  INTRODUCTION 

2.1  Background  on  moment  tensor  source  type 

The  moment  tensor  is  an  expression  of  the  source  for  small  to  moderate  seismic 
events.  Mathematically,  it  is  a  3  x  3  symmetric  matrix,  and  the  space  of  moment  tensors 
is  therefore  six-dimensional.  Traditionally,  moment  tensors  have  often  been  constrained 
to  be  deviatoric,  meaning  that  their  trace  should  be  zero,  in  which  case  only  a  subspace  of 
the  six-dimensional  space  is  used.  In  the  past  two  decades,  however,  several  authors  have 
argued  against  such  constraints,  especially  when  the  seismic  events  under  consideration 
are  not  true  tectonic  earthquakes  but  rather  are  associated  with  volcanic  or  geothermal 
activity,  glaciers,  mine  collapses,  hydraulic  fracturing,  or  explosions  (Ford  et  al.,  2009, 
2010);  see  review  in  Tape  and  Tape  (2013,  Figures  14  and  S14). 

The  emphasis  within  our  proposed  study  is  on  the  uncertainties  of  moment  tensors, 
notably  those  associated  with  the  source  type.  The  source  type  refers  to  two  coordinates 
that  characterize  the  isotropic  and  CLVD  components  within  a  moment  tensor.  These 
components  are  provided  by  the  lune  latitude  (5)  and  lune  longitude  (y),  which  are  simply 
angular  measures  in  eigenvalue  space.  The  ‘fundamental  lune’  is  a  sector  of  this  3D 
eigenvalue  space  that  covers  all  possible  moment  tensor  source  types  (cf.  Riedesel  and 
Jordan,  1989).  The  center  of  the  lune  is  a  double  couple  (1,  0,-1),  while  the  top  and 
bottom  of  the  lune  are  isotropic  points  (Figure  1).  The  lune  can  be  partitioned  into  four 
regions  by  three  arcs  at  fa  =  0.  'ki  =  0,  and  A]  =  0.  The  distance  on  the  lune  is  the  natural 
measure  of  distance  between  moment  tensors.  For  example,  for  two  moment  tensors  with 
the  same  basis  U,  the  angle  between  lune  points  is  the  same  as  the  angle  between  the 
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moment  tensors  in  matrix  space  (Tape  and  Tape,  2013).  The  ‘amount’  of  double  couple 
or  CLVD  or  isotropic  part  within  a  moment  tensor  is  calculated  as  the  angular  distance 
from  a  moment  tensor  to  one  of  these  points  on  the  lune. 

The  lune — by  which  we  mean  the  eigenvalue  space  for  moment  tensors — provides 
clarity  to  some  convoluted  areas  of  moment  tensor  literature.  For  example,  there  are  three 
different  published  formulas  for  the  ‘percent’  of  ISO,  DC,  and  CLVD  from  a  given 
moment  tensor;  furthermore,  two — not  three — parameters  should  be  needed  to 
characterize  moment  tensor  source  type.  And  there  are  at  least  four  published  equations 
for  obtaining  the  seismic  moment  from  a  given  moment  tensor;  see  comparison  in  Tape 
and  Tape  (2012b,  Figure  13).  Thus  it  is  challenging  to  compare  results  across  studies 
when  so  many  different  methods  are  used  to  represent  the  solutions.  Finally,  it  turns  out 
that  alternate  representations  for  source  type — notably  the  T-k  plot  of  Hudson  et  al. 
(1989) — turn  out  to  be  underlain  by  unphysical  assumptions  about  uniformity  of  source 
types  (Tape  and  Tape,  2012b).  The  lune  provides  a  natural  framework  that  leads  to 
simpler  expressions  for  characterizing  moment  tensors. 

Three  relevant  topics  from  Tape  and  Tape  (2013)  are  (1)  the  calculation  of  volume 
change,  (2)  the  lack  of  justification  for  the  deviatoric  constraint,  and  (3)  the  interpretation 
of  multiple  events.  Tape  and  Tape  (2013)  called  into  question  the  ubiquitous  (in  the 
published  literature)  interchangeability  of  ‘isotropic  component’  and  ‘volumetric 
component’.  In  fact,  the  volume  change  component  varies  with  longitude — not  latitude — 
on  the  lune.  We  also  questioned  the  continued  use  of  the  deviatoric  constraint  within 
operational  moment  tensor  codes:  this  constraint  is  only  needed  because  it  is  difficult  to 
constrain  uncertainties.  (In  fact,  from  a  theoretical  standpoint,  the  constant-n  constraint 
would  be  more  appropriate  than  the  zero-trace  constraint  (Minson  et  al.,  2007;  Tape  and 
Tape,  2012a).)  We  also  introduced  a  generic  ‘two-source  model’  that  could  produce  a 
wide  spread  of  full  moment  tensors  on  the  lune  that  is  seen  in  the  literature.  But  in 
general  we  caution  against  invoking  multiple  events  or  using  decompositions,  as  the 
classical  model  for  moment  tensors — an  oblique  opening  (or  closing)  crack  within  a 
medium  with  Poisson  parameter  u — provides  the  simplest  interpretation  for  full  moment 
tensors  (Tape  and  Tape,  2013).  As  noted  in  Tape  and  Tape  (2013),  we  must  have 
uncertainty  estimates  for  published  moment  tensors  in  order  to  test  physical  models  for 
moment  tensors. 

2.2  Theoretical  development  and  application  for  moment  tensor  uncertainties 

In  the  past  decade,  efforts  have  been  made  to  characterize  uncertainties  of  moment 
tensors.  Some  studies  have  estimated  the  uncertainties  by  repeatedly  solving  for  Mo  using 
different  subsets  of  data  (boot- strapping  or  related  Monte  Carlo  approaches).  Other 
studies  have  provided  a  view  of  how  the  misfit  function  varied  over  the  full  space  of 
moment  tensor  source  types  (Ford  et  al.,  2010;  Alvizuri  and  Tape,  2016).  To  our 
knowledge,  Stahler  and  Sigloch  (2014)  is  the  only  study  to  have  formally  accounted  for 
uncertainties,  by  generating  moment  tensor  samples  of  a  posterior  probability  density. 
This  was  made  possible  by  using  a  uniform  distribution  of  moment  tensors;  we  use  an 
alternative  uniform  distribution  in  our  study  Tape  and  Tape  (2015  a). 

The  goal  of  this  study  is  to  estimate  moment  tensor  uncertainties  for  three  sets  of 
previously  studied  events:  21  Alaska  earthquakes,  63  volcanic  events,  and  29  events  from 
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the  Nevada  Test  Site.  Following  Tape  and  Tape  (2016),  we  represent  moment  tensor 
uncertainties  in  the  form  of  a  confidence  curve  that  expresses  the  concentration  of 
probability  in  moment  tensor  space  near  Mo.  We  build  upon  three  recent  studies:  Tape 
and  Tape  (2016),  which  established  the  theory  for  confidence  curves;  Silwal  and  Tape 
(2016),  which  estimated  confidence  curves  for  double  couple  moment  tensors;  and 
Alvizuri  and  Tape  (2016),  which  estimated  a  catalog  of  full  moment  tensors  for  volcanic 
events.  In  our  examination  of  earthquakes  and  nuclear  explosions  from  the  Nevada  Test 
Site,  we  use  the  exact  same  events  as  in  Ford  et  el.  (2009).  In  comparison  with  the 
approach  of  Ford  et  el.  (2009),  we  use  more  stations,  body  waves  (in  addition  to  surface 
waves),  a  broader  bandpass  for  surface  waves  (T  =  7-50  s),  and  an  alternative 
characterization  of  moment  tensor  uncertainties. 

2.3  Summary  of  project  publications 

Our  project  resulted  in  eight  publications,  including  three  publications  in  peer-reviewed 
journals: 

•  Journal  publications:  Tape  and  Tape  (2015a,  2016,  2017) 

•  Journal  publication  in  preparation:  Alvizuri  et  al.  (2017,  in  prep.) 

•  Abstracts:  Tape  and  Tape  (2015b);  Alvizuri  et  al.  (2016,  2017);  Tape  (2017a) 

•  Technical  reports:  Tape  (2017b) 

We  expect  to  submit  a  fourth  publication  (Alvizurti  et  al.,  2017,  in  prep.)  by  December 
2017.  Our  final  report  focuses  on  the  unsubmitted  manuscript,  since  full  details  of  the 
other  work  can  be  found  within  the  three  published  papers. 

3.  TECHNICAL  APPROACH 

3.1  Data 

We  examine  events  from  three  previous  studies:  21  earthquakes  in  southern  Alaska 
(tab.  4,  Silwal  and  Tape,  2016),  63  events  at  Uturuncu  volcano  in  Bolivia  (tab.  4,  Alvizuri 
and  Tape,  2016),  and  29  earthquakes  and  nuclear  tests  at  the  Nevada  Test  Site  (tab.  1  of 
Ford  et  al.,  2009).  Figure  1  shows  the  three  regions  and  the  stations  used  within  the 
moment  tensor  inversions. 

For  the  Alaska  and  Uturuncu  data  sets,  all  waveforms  are  openly  available  from  the 
IRIS  Data  Management  Center.  For  the  Nevada  Test  Site  events,  the  waveforms  are 
openly  available  from  IRIS  and  the  Northern  California  Data  Center,  while  other 
waveforms  are  from  Walter  et  al.  (2006). 

All  waveforms  were  downloaded  and  processed  using  ObsPy,  a  python-based 
package  for  seismology.  The  processing  steps  for  each  event  were:  (1)  obtain  three- 
component  waveforms  and  metadata  from  the  IRIS  (Incorporated  Research  Institutions 
for  Seismology)  Data  Management  Center,  the  Northern  California  Earthquake  Data 
Center  (NCEDC),  or  Walter  et  al.  (2006);  the  time  interval  is  100~s  before  the  origin 
time  to  600~s  after  the  origin  time,  (2)  remove  instrument  response  using  an  acausal 
Butterworth  filter,  (3)  using  the  source-station  azimuth  and  the  sensor  orientation  angle, 
rotate  horizontal  components  to  radial  and  transverse  directions.  Additional  processing 
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steps,  such  as  cutting  time  windows  and  additional  bandpass  filtering,  were  applied 
during  the  moment  tensor  inversions. 

We  use  high-frequency  P  waves  for  the  Uturuncu  and  NTS  events.  Within  the  misfit 
function,  described  next,  we  make  waveform  measurements  between  observed  and 
synthetic  P  waves.  These  measurements  are  more  reliable  when  we  align  the  P  waves  on 
the  observed  P  arrival  time. 

3.2  Methods 

Seismic  moment  tensor  inversion  requires  specifying  a  misfit  function  between 
observed  seismograms  and  synthetic  (or  'modeled')  seismograms.  Synthetic  seismograms 
are  calculated  using  a  model  of  the  source— in  our  case,  a  moment  tensor— and  a  model 
of  Earth's  structure— in  our  case,  a  layered  model.  We  calculate  synthetic  seismograms 
for  a  range  of  source  depths  and  source-station  distances  using  the  frequency- 
wavenumber  method.  The  source  origin  times  and  hypocenters  are  assumed  to  be  fixed. 

The  misfit  function  we  use,  from  Zhu  and  Helmberger  (1996),  measures  integrated 
differences  between  observed  and  synthetic  velocity  seismograms  within  five  time 
windows:  P  wave  on  vertical  and  radial  components,  Rayleigh  wave  on  vertical  and 
radial  components,  and  Love  wave  on  transverse  component.  Prior  to  calculating  the 
waveform  difference,  we  allow  the  synthetic  seismograms  to  be  shifted  in  time  to 
maximize  the  cross-correlation  with  observed  seismograms.  For  each  station,  three  time 
shifts  are  allowed,  for  the  P,  Rayleigh,  and  Love  waves.  These  different  time  shifts  are 
allowed  because  these  waves  are  sensitive  to  different  parts  of  Earth  structure.  By 
allowing  for  different  time  shifts,  we  are  recognizing  that  our  layered  Earth  structural 
models  cannot  adequately  predict  the  traveltimes  for  seismic  waves.  Additional 
discussion  of  time  shifts  can  be  found  in  Silwal  and  Tape  (2016). 

For  our  moment  tensor  inversions,  we  use  the  ‘cut-and-paste’  (CAP)  code  of  Zhu 
and  Helmberger  (1996)  and  Zhu  and  Ben-Zion  (2013),  with  some  modifications.  For  the 
misfit  function  we  use  an  LI  norm  Silwal  and  Tape  (2016),  and  we  incorporate  the 
number  of  misfitting  polarities  into  the  waveform-based  misfit  function.  We  also  use  an 
efficient  grid  search  over  moment  tensor  space,  described  below. 

3.3  Moment  tensor  source  type 

The  source  type  of  a  seismic  moment  tensor  is  a  represented  by  two  parameters  that 
are  derived  directly  from  the  three  eigenvalues  of  the  3x3  moment  tensor.  The  norm  of 
the  eigenvalue  triple  provides  the  magnitude  of  the  moment  tensor,  leaving  two  free 
parameters  to  define  the  source  type.  In  the  same  year,  Hudson  et  al.  (1989)  introduced 
the  T-k  source  type  plot,  while  Riedesel  and  Jordan  (1989)  introduced  a  representation  of 
source  types  on  a  sphere.  The  motivation  for  the  2D  T-k  plot  was  to  achieve  uniform 
distribution  of  moment  tensors — however  uniformity  was  assumed  to  apply  to  a  cube  in 
eigenvalue  space,  an  assumption  with  no  reasonable  justification.  Chapman  and  Leaney 
(2012)  dismissed  the  fundamental  assumption  of  the  T-k  representation,  and  they 
advocated  using  the  eigensphere,  following  Riedesel  and  Jordan  (1989).  Tape  and  Tape 
provided  a  geometrical  framework  for  using  the  eigensphere  for  representing  source  type, 
with  the  two  parameters  being  lune  longitude — quantifying  the  CLVD  component — and 
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lune  latitude — quantifying  the  isometric  component.  Tape  and  Tape  (2012b)  compared 
the  source  type  representations  of  T-k  and  the  lune,  showing  the  problems  associated  with 
the  cube  assumption  underlying  the  T-k  representation. 

As  pointed  out  by  several  authors,  including  Hudson  et  al.  (1989),  Ford  et  al.  (2010), 
and  Chapman  and  Leaney  (2012),  neither  the  source  type  plots  of  T-k  nor  the  lune 
provide  a  uniform  distribution  of  moment  tensors.  Toward  the  boundary  of  the  source 
type  plots,  the  eigenvalues  become  degenerate  and  the  moment  tensors  are  similar  to  each 
other,  irrespective  of  orientation.  Tape  and  Tape  (2015)  introduced  a  parameterization  of 
source  type  that  provides  uniform  moment  tensors:  “The  parametrization  is  uniform,  in 
the  sense  that  equal  volumes  in  the  coordinate  domain  of  the  parametrization  correspond 
to  equal  volumes  of  moment  tensors.”  (Tape  and  Tape,  2015,  used  u  and  v,  but  in  our 
work  we  have  used  w  instead  of  u;  u  is  similar  to  colatitude,  while  w  is  similar  to 
latitude.)  Tape  and  Tape  (2015)  discuss  the  vw  rectangle  in  comparison  with  the  lune: 

“The  homogeneous  probability  for  moment  tensors  is  the  probability  in  which  equal 
volumes  of  moment  tensors  are  equally  likely.  If  we  are  committed  to  the  homogeneous 
probability  as  the  a  priori  probability  for  moment  tensors,  then  the  lune  is  not  the  right 
place  to  display  eigenvalue  triples  0,  if  the  display  is  intended  to  reveal  non-uniformities. 
(For  example,  a  tendency  for  the  0s  to  favour  the  center  of  the  lune  does  not  mean  that  we 
are  dealing  with  a  double  couple.)  Rather  the  eigenvalue  triples  should  be  displayed  in 
the  uv  coordinate  domain  Q,  since  the  a  priori  probability  density  for  eigenvalue  triples  in 
Q  is  constant.  When  not  specifically  trying  to  convey  or  detect  non-uniformities, 
however,  we  recommend  the  lune,  since,  after  all,  the  lune  is  exactly  the  set  of  ordered 
and  normalized  eigenvalue  triples.” 

We  see  merits  to  both  the  lune  and  the  vw  plot.  We  see  no  justification  for  the  underlying 
assumption  associated  with  the  T-k  plot;  nevertheless,  the  plot  still  serves  the  purpose  of 
conveying  the  moment  tensor  source  type. 

4.  RESULTS 

4.1  Nevada  earthquakes 

The  Little  Skull  Mountain  earthquake  of  1992-06-29  occurred  inside  the  Nevada 
Test  Site  (Figure  2)  and  is  the  largest  earthquake  analyzed  by  Ford  et  al.  (2009).  It  is 
notable  for  its  size  (Mw  >  5)  and  its  occurrence  near  Yucca  Mountain,  which,  at  the  time, 
was  under  consideration  as  storage  site  for  nuclear  waste.  Table  1  of  Lohman  et  al. 
(2002)  summarizes  several  studies  of  the  source  parameters  of  the  Little  Skull  Mountain 
earthquake. 

Our  moment  tensor  inversions  for  the  Little  Skull  Mountain  earthquake  used  15 
broadband  stations.  For  our  first  inversion,  we  assume  a  double  couple  moment  tensor 
and  include  body  waves  (filtered  1-4  s)  and  surface  waves  (filtered  10-50  s).  We  first 
perform  a  coarse  search  over  magnitude,  depth,  strike,  dip,  and  rake,  to  determine  the 
best-fitting  depth  and  magnitude.  The  depth  search  over  full  moment  tensor  space 
provides  similar  results,  with  a  best-fitting  depth  of  10  km.  We  use  10  km  for  the 
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estimated  depth,  based  on  consensus  from  previous  studies,  several  of  which  are  more 
sensitive  to  depth  estimation  (Lohman  et  al.,  2002). 

With  the  depth  and  magnitude  fixed,  we  then  perform  a  fine  grid  search  over  the 
space  of  strike,  dip,  and  rake.  Synthetic  seismograms  from  the  ID  model  are  able  to  fit 
the  main  surface  wave  arrivals  on  the  data,  but  they  are  not  able  to  fit  the  complexities 
caused  by  3D  structure,  such  as  for  paths  across  the  San  Joaquin  Valley  and  Sierra 
Nevada.  These  effects  are  diminished  when  filtering  seismograms  with  a  larger  minimum 
period  (such  as  20  s  instead  of  10  s). 

Our  misfit  analysis  of  double  couple  moment  tensors  was  presented  in  Silwal  and 
Tape  (2016)  for  Alaska  earthquakes.  Application  to  the  Little  Skull  Mountain  earthquake 
provides  a  sampling  of  20  beachballs  from  the  posterior  distribution;  we  see  that  a 
normal-fault  mechanism  is  preferred,  but  that  the  data  allow  for  some  notable  deviations. 
An  optimal  confidence  curve  is  one  that  has  a  steep  slope  at  the  origin  and  then  rises  all 
the  way  to  P(V)  ~  1,  such  that  P(V)  makes  the  shape  of  a  L  For  the  Little  Skull  Mountain 
earthquake,  the  confidence  curve  has  a  steep  slope  at  the  origin  but  then  reveals  a  more 
gradual  rise,  indicating  that  the  probability  for  the  solution  is  encountered  gradually  as 
one  moves  away  from  the  global  minimum  moment  tensor  solution. 

The  full  moment  tensor  solution  for  Little  Skull  Mountain  reveals  a  minor 
improvement  of  waveform  fits  (VR  30.6  vs  28.2).  Our  initial  misfit  analysis  of  full 
moment  tensors  was  presented  in  Alvizuri  and  Tape  (2016)  for  events  at  Uturuncu 
volcano.  Application  to  the  Little  Skull  Mountain  earthquake  reveals  a  best-fitting 
solution  that  is  close  to  the  deviatoric  arc,  with  lune  longitude  y  =  14°  and  lune  latitude  3°. 

Within  the  moment  tensor  inversions,  time  shifts  are  applied  to  the  synthetic 
waveforms  in  an  attempt  to  account  for  the  approximateness  of  our  layered  velocity 
model.  The  time  shifts  provide  a  perspective  on  potential  cycle-skipping  of  waveform 
fits,  since  stations  in  the  same  azimuth  should  have  time  shifts  with  the  same  sign. 

The  time  shifts  exhibit  systematic  variations  with  azimuth  that  we  attribute  to 
differences  in  3D  structure.  The  sign  discrepancy  between  time  shifts  for  Love  and 
Rayleigh  waves  suggests  the  presence  of  crustal  anisotropy  or  that  the  ID  model  needs  to 
be  adjusted  in  the  shallowest  layer,  where  the  Love  and  Rayleigh  waves  have  differing 
sensitivities.  The  plots  also  reveal  stations  that  are  nodal  to  the  mechanism.  At  these 
stations,  the  P  and  Rayleigh  waves  can  be  low  amplitude  and  possibly  distorted  by  3D 
structure  effects. 

4.2  Nevada  nuclear  tests 

Following  Ford  et  al.  (2009),  we  use  the  1991-09-14  HOYA  nuclear  test  to  illustrate 
the  waveform  fits  and  misfit  analysis.  For  HOYA  we  use  14  stations  and  filter  surface 
waves  at  10-50  s  and  body  waves  at  high  frequencies  (0.8-10  Hz).  For  the  LLNL  stations 
we  use  the  ‘HF’  waveforms  for  the  body  waves  and  ‘LF’  waveforms  for  the  surface 
waves.  (These  waveforms  are  recorded  by  different  sensors.)  For  the  body  waves,  we 
align  the  synthetic  P  waves  on  the  observed  P  onset  times  and  then  allow  a  time  shift  of 
only  0.5  s. 

The  waveform  fits  and  misfit  analysis  for  HOYA  are  shown  in  Figures  3  and  4.  The 
best-fitting  moment  tensor  is  at  (5°,  76  °)  on  the  lune,  within  the  isotropic  region  defined 
by  all-positive  eigenvalues.  The  region  of  lowest  fit  is  tightly  constrained  around  the  best- 
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fitting  solution,  but  there  is  also  a  weak  local  minimum  in  that  negative  isotropic  region 
(Figure  4)  arising  from  cycle- skipping. 

We  calculate  moment  tensor  solutions  for  all  32  events  from  Ford  et  al.  (2009).  Lune 
plots  and  waveform  fits  for  four  nuclear  explosions  are  shown  in  Figures  5  and  6.  The 
events  of  Ford  et  al.  (2009)  are  represented  as  beachballs  on  the  lune  in  Figure  9a.  We 
also  plot  the  beachballs  on  the  vw  rectangle,  which  has  an  advantage  over  the  lune  in  that 
distances  on  the  plot  better  reflect  distances  between  moment  tensors.  As  an  example,  the 
17  explosions  from  Ford  et  al.  (2009)  are  spread  out  on  the  lune  but  are  close  together  on 
the  rectangle,  indicating  that  the  moment  tensors  are  quantifiably  very  similar  to  each 
other.  Our  moment  tensor  solutions  for  the  32  events  are  shown  in  Figure  9b. 

4.3  Uturuncu  volcanic  events 

We  reproduce  the  results  from  Alvizuri  and  Tape  (2016),  but  using  a  modified  misfit 
function.  Previously  we  excluded  any  moment  tensor  that  had  any  disagreement  between 
its  predicted  first-motion  polarities  and  the  observed  first-motion  polarities.  The  modified 
misfit  function  balances  the  polarity  misfit  with  the  waveform  misfit.  This  means  that  the 
misfit  function  is  plotted  on  the  entire  lune,  rather  than  on  a  subregion,  as  in  Alvizuri  and 
Tape  (2016). 

Figure  9c  shows  the  full  moment  tensor  catalog  for  63  small  (Mw  <  3)  events  at 
Uturuncu  volcano.  Most  events  in  the  catalog  have  a  positive  isotropic  component  (8  > 
0°).  We  show  the  beachballs  on  a  plot  of  the  vw  rectangle,  which  is  more  appropriate  than 
the  lune  when  the  intent  is  to  represent  distances  in  moment  tensor  space.  All-outward 
beachballs  (with  all  positive  eigenvalues)  are  quite  similar  to  each  other;  as  we  might 
expect,  they  occupy  a  much  larger  region  on  the  lune  than  they  do  on  the  vw  rectangle. 

4.4  Alaska  earthquakes 

We  apply  the  analysis  of  Alvizuri  and  Tape  (2016)  to  the  21  earthquakes  of  Silwal 
and  Tape  (2016).  No  first-motion  polarities  were  used,  so  the  misfit  function  uses  only 
waveform  differences.  Figure  9d  shows  the  full  moment  tensor  solutions  on  the  lune  and 
on  the  vw  rectangle.  While  some  earthquakes  exhibit  a  possibly  spurious,  negative 
isotropic  component,  most  events  are  quite  close  to  the  double  couple  at  the  center  of  the 
lune. 

5.  DISCUSSION 

Earthquakes  are  useful  for  calibrating  the  results  from  non-earthquake  events,  such 
as  nuclear  explosions.  Importantly,  the  full  moment  tensor  inversions  for  earthquakes 
need  to  be  performed  with  conditions  that  are  as  close  to  possible  to  those  of  the  non- 
earthquake  events,  notably:  same  stations,  same  time  period  (meaning  same  station 
metadata),  similar  epicentral  region  (similar  paths).  This  was  the  motivation  behind  the 
event  selection  in  Ford  et  al.  (2009). 

Earthquake  mechanisms  estimated  as  full  moment  tensors  tend  to  cluster  near  the 
double  couple  region  of  source  type  space,  as  we  might  hope.  However,  the  spread  of 
best-fitting  moment  tensors  can  deviate  significantly  from  the  ideal  double  couple,  as 
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shown  in  Boyd  et  al.  (2015)  for  a  set  of  828  events  in  the  Berkeley  Seismological 
Laboratory  catalog.  The  Berkeley  moment  tensor  solutions  typically  use  <8  stations  and 
are  based  on  low-pass-filtered  surface  waves.  By  comparison,  our  full  moment  tensor 
analysis  of  21  Alaska  earthquakes  uses  dozens  of  stations  and  also  includes  body  waves. 
The  moment  tensors  are  more  tightly  clustered  around  the  double  couple  at  the  center  of 
the  lune  (Figure  9d).  This  suggests  that  improved  data  coverage  and  use  of  body  waves 
can  diminish  the  spread  in  source  types  that  are  observed  for  sets  of  earthquakes. 

For  the  moment  tensor  inversions  in  this  study  we  have  included  waveform  fits  for 
body  waves.  To  examine  the  influence  of  body  waves,  we  performed  separate  inversions 
with  body  waves  only,  surface  waves  only,  and  both  combined.  Results  for  the  Little 
Skull  Mountain  earthquake  and  the  HOYA  nuclear  test  are  shown  in  Figure  8.  For  the 
Little  Skull  Mountain  earthquake  we  see  that  all  three  mechanisms  are  notably  different, 
though  they  all  share  a  white  central  region  and  solid  regions  to  the  northwest  and 
southeast.  The  body+surf  lune  plot  (Figure  8c)  appears  to  be  a  blend  of  the  lune  plots  for 
the  body-only  (Figure  8a)  and  surface-only  (Figure  8b)  cases.  Note  that  the  magnitude 
estimate  from  body  waves  only  is  lower  than  the  overall  estimate  (Mw  4.9  vs  Mw  5.2). 

For  HOYA,  the  lune  plot  for  body+surf  (Figure  8f)  reveals  multiple  good- fitting 
regions,  some  of  which  are  present  in  the  body-only  (Figure  8d)  and  surface-only  (Figure 
8e)  lune  plots.  Although  first-motion  polarities  are  not  used  within  the  inversion,  each 
mechanism  predicts  outward  first-motions  at  all  stations. 

It  is  possible  that  the  mechanisms  obtained  from  high-frequency  body  waves  could 
differ  from  the  mechanisms  obtained  from  long-period  surface  waves.  This  could  arise 
from  time-dependent  complexity  in  the  source  process,  whereby  the  initial  slip  is 
responsible  for  the  body  waves  and  later,  longer-duration  slip  is  responsible  for 
generating  surface  waves.  Therefore  we  must  be  careful  in  interpreting  differences  in 
Figure  8,  since  some  could  be  real,  rather  than  caused  by  errors  in  modeling  or  data. 

Our  magnitude  estimates  for  the  Little  Skull  Mountain  earthquake  and  for  HOYA  are 
lower  than  those  in  the  published  literature.  Our  Mw  5.2  estimate  for  Little  Skull 
Mountain  compares  with  Mw  5.6  from  Lohman  et  al.  (2002)  and  Ford  et  al.  (2009)  and 
Mw  5.7  from  the  GCMT  catalog.  Further  work,  including  testing  of  different  velocity 
models  and  attenuation  models,  is  needed  to  examine  this  discrepancy. 

The  application  of  time  shifts  to  synthetic  seismograms  in  moment  tensor  inversions 
is  a  widely  used  practice  that  can  be  both  subtle  and  dangerous.  As  an  extreme  example, 
the  moment  tensor  -Mo  will  exhibit  the  same  relative  variations  in  amplitudes  as  the  best¬ 
fitting  moment  tensor  Mo.  While  the  inclusion  of  first-motion  polarity  measurements  can 
distinguish  between  these  extreme  cases,  there  are  other  cases  that  are  difficult.  With  only 
ten  or  so  waveforms,  as  in  the  case  of  the  NTS  events,  a  single  misaligned  waveform  can 
influence  the  solution.  One  challenge  is  that  the  prospects  of  cycle  skipping — 
misalignment  of  the  synthetic  wave  with  the  observed  wave — increase  as  we  decrease  the 
period  content. 

The  best  approach  to  eliminating  time  shifts  would  be  to  use  a  3D  velocity  model 
that  accurately  predicts  the  arrival  times  of  all  waves  used  in  the  moment  tensor 
inversion.  For  us,  this  would  require  accurate  prediction  of  the  seismic  wavefield  for 
periods  >1  s.  Such  a  regional  model  does  not  exist  for  any  active  tectonic  setting, 
including  the  western  United  States.  Cycle  skipping  of  surface  waves  is  particularly 
problematic,  so  one  could  have  the  more  realistic  goal  of  eliminating  time  shifts  for 
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periods  >5  s.  With  dense  data  arrays  and  significant  computational  modeling,  it  is 
possible  to  fit  the  seismic  wavefield  at  these  periods.  To  integrate  such  a  model  into  a 
moment  tensor  inversion  would  require  replacing  the  inexpensive,  efficient  calculation  of 
ID  synthetic  seismograms  with  expensive  calculation  of  3D  synthetic  seismograms  from 
wavefield  simulations. 

6.  CONCLUSIONS 

We  present  full  moment  tensor  solutions,  with  uncertainties,  for  three  sets  of  events: 
21  earthquakes  in  southern  Alaska,  63  events  at  Uturuncu  volcano,  Bolivia,  and  29  events 
at  the  Nevada  Test  Site.  We  characterize  each  event  as  a  point-source  moment  tensor 
with  a  fixed  epicenter  and  origin  time.  The  source  time  function  is  a  trapezoidal  function 
whose  width  scales  with  magnitude.  We  perform  a  separate  grid  search  over  depth  and 
magnitude,  and  our  emphasis  on  uncertainty  estimation  is  with  the  5D  space  of 
normalized  moment  tensors. 

Our  main  findings  are  as  follows: 

1.  The  misfit  function  in  moment  tensor  space  exhibits  local  minima  and  other 
complexities  (e.g.,  Figure  4).  These  reveal  the  importance  of  evaluating  the  misfit 
function  over  the  full  space  of  moment  tensors,  as  opposed  to  linearized 
inversions,  which  may  lead  into  an  incorrect  local  minimum. 

2.  The  complexities  of  the  misfit  function  can  be  distilled  into  a  confidence  curve 
P(V)  that  adjusts  the  posterior  probability  for  the  homogeneous  probability  of 
uniform  moment  tensors  (Tape  and  Tape,  2016;  Silwal  and  Tape,  2016).  The  area 
under  the  confidence  curve  is  the  confidence  parameter  Pav.  Even  in  cases  where 
moment  tensor  space  is  uniformly  parameterized  (e.g.,  Stahler  and  Sigloch,  2014), 
there  is  still  a  benefit  to  the  confidence  curve  in  that  it  can  account  for  the 
‘imprinted’  influence  of  the  homogeneous  probability  for  moment  tensors. 

3.  We  present  a  high-quality  catalog  of  full  moment  tensors  from  the  Nevada  Test 
Site,  building  upon  the  efforts  of  Ford  et  al.  (2009).  In  comparison  with  Ford  et  al. 
(2009),  we  include  more  stations,  we  use  a  lower  minimum  period  (10  s)  for  all 
events,  and  we  include  (regional)  body  waves  in  all  inversions.  The  inclusion  of 
body  waves  improves  the  coverage  of  the  source  hemisphere,  since  the  body- 
wave  take-off  angles  tend  to  be  steep  (and  downward),  whereas  surface-wave 
paths  are  near  horizontal.  Further  efforts  to  include  first-motion  polarity 
measurements  from  short-period  stations  could  result  in  moment  tensor  solutions 
with  higher  confidence  than  those  presented  here. 

Additional  work  is  needed  to  improve  our  understanding  of  the  theory  and  application  of 
moment  tensor  confidence  curves.  Most  choices  are  made  within  the  misfit  function, 
which  includes  time  windowing,  weighting,  choice  of  norm,  and  time  shifting  of 
synthetic  seismograms.  For  example,  a  more  complete  representation  of  data  covariance 
terms  could  eliminate  the  need  for  a  scale  factor.  The  use  of  3D  synthetic  seismograms 
could  lessen  the  need  for  time  shifts,  if  the  3D  velocity  model  is  accurate.  Some  questions 
can  be  adequately  addressed  without  real  data  within  synthetic  inversions. 
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Our  primary  goal  has  been  to  characterize  uncertainties  of  moment  tensors.  This  provides 
a  starting  point  for  interpreting  the  results  in  the  context  of  physical  models,  such  as  those 
discussed  in  Tape  and  Tape  (2013)  or  Patton  and  Taylor  (201 1).  We  can  interpret  any  full 
moment  tensor  as  a  one-process  model  of  an  oblique  opening  crack  (Minson  et  al.,  2007) 
or  as  a  two-process  model  of  a  tensional  or  compressional  crack  combined  with  shear  slip 
within  the  same  crack  plane  (Aki  and  Richards,  1980).  More  complex  source  processes 
have  been  proposed,  but  they  cannot  be  distinguished  or  evaluated  using  the  typical  set  of 
seismic  waveforms. 
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Figure  1.  Fundamental  lune  for  moment  tensor  source  types.  Four  regimes  for  moment 
tensors  on  the  fundamental  lune  L.  Above  and  to  the  right  of  the  A2  =  0  (blue  curve), 
moment  tensor  beachballs  have  red  bands  and  white  caps  (A2  >  0  >  A3).  Below  and  to  the 
left  of  it  they  have  white  bands  and  red  caps  (Ai  >  0  >  A2).  Above  the  arc  A3  =  0  (red), 
beachballs  are  all  red  (A3  >  0).  Below  the  Ai  =  0  (white),  they  are  all  white  (0  >  Ai). 
(from  Tape  and  Tape,  2013,  Figure  SI) 
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Figure  2.  Locator  maps  for  study  regions.  Three  study  regions  of  the  Nevada  Test  Site, 
Alaska,  and  Bolivia,  (a)  Global  perspective  showing  our  three  study  regions.  ( b )  Western 
United  States,  showing  Nevada  Test  site  (outlined),  (c)  Zoom-in  on  Nevada  Test  Site. 
Nuclear  explosions  (17)  are  red,  earthquakes  (12)  are  blue,  and  collapses  (3)  are  yellow, 
(d)  21  earthquakes  in  southern  Alaska,  (e)  63  events  at  Uturuncu  volcano,  Bolivia. 


Approved  for  public  release;  distribution  is  unlimited. 
12 


Figure  3.  Waveform  fits  for  HOYA  nuclear  test. 
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Moment  tensor  solution  and  subset  of  waveform  comparisons  for  the  HOYA  nuclear 
explosion.  Each  column  is  a  different  section  of  the  three-component  waveform:  PV  = 
vertical  component  P  wave,  PR  =  radial  component  P  wave,  SurfV  =  vertical  component 
Rayleigh  wave,  SurfR  =  radial  component  Rayleigh  wave,  SurfT  =  transverse  component 
Love  wave.  The  stations  are  ordered  by  increasing  epicentral  distance  from  the  top  row. 
The  observed  waveforms  are  plotted  in  black,  the  synthetic  waveforms  are  plotted  in  red. 
The  body  waves  are  filtered  3-10  Hz,  and  the  surface  waves  are  filtered  7-33  s.  The 
numbers  below  each  station  name  are  the  station  epicentral  distance  (top)  and  station 
azimuth  (bottom).  The  four  numbers  below  each  pair  of  waveforms  are,  from  top  to 
bottom,  (1)  the  cross-correlation  time  shift  AT  =  Tots  ~  Tyn  required  for  matching  the 
synthetics  s(t)  with  the  data  u(t)  (a  positive  time-shift  means  that  the  synthetics  arrive 
earlier  than  the  data);  (2)  the  maximum  cross-correlation  percentage  between  u(t)  and 
s(t  -  AT);  (3)  the  percentage  of  the  total  misfit;  and  (4)  the  amplitude  ratio  In(Aobs/ASyn) 
in  each  time  window. 
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Event  HOYA,  M  4.70 
Lon -116.4290,  Lat  37.2260 
Dep  0.7  km  (inversion  1  km) 
Best  source  type  y  5°,  6  76° 


(1,1  ,-2) 

CLVD 


Figure  4.  Uncertainty  summary  for  HOYA  nuclear  test. 
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Waveforms  for  the  best-fitting  moment  tensor  are  shown  in  Figure  3.  (a)  Map  of  source 
location  (red  star)  and  stations  used  in  the  inversion  for  this  event,  (b)  Contour  plot  of  the 
polarity  misfit  on  the  lune,  if polarities  are  used ,  as  in  the  case  of  the  Uturuncu  events 
(Alvizuri  and  Tape,  2016).  (c)  Probability  density  for  moment  tensor  source  type,  plotted  in  vw 
space  (Tape  and  Tape,  2015).  A  green  circle  indicates  the  location  of  the  maximum,  (d)  Contour 
plot  of  the  variance  reduction  VR(A)  .  The  variance  reduction  VR(A)  at  a  point  A 
is  the  maximum  variance  reduction  VR(M)  for  moment  tensors  M  that  have  source  type  A 
and  that  have  correct  polarities.  Large  values  (blue)  ofVR  represent  better  fit  between 
observed  and  synthetic  waveforms.  Of  the  beachballs  M(A),  our  desired  solution  Mo  (green 
box)  is  the  one  with  largest  VR.  The  gray  arcs  on  the  lune  are  the  great  circle  arcs  X  /  =  0, 

A 2  =  0,  and  2.3  =  0  (white,  cyan,  and  red  in  Figure  1).  Selected  eigenvalue  triples  (black 
dots)  on  the  boundary  of  the  lune  are  indicated,  with  the  understanding  that  the  triples 
need  to  be  normalized.  The  positive  isotropic  source  (1,  1,  1)  is  at  the  top,  the  negative 
isotropic  source  (—1—1—1)  is  at  the  bottom,  and  the  double  couple  (1,  0,-1),  not  shown, 
would  be  at  the  center  of  the  lune.  (e)  Uniform  and  posterior  probability  curves  with 
respect  to  co,  the  angular  distance  from  Mo  to  a  moment  tensor,  (f)  Confidence  curve  P(V) 
showing  the  distribution  of  probability  as  a  function  of  fractional  volume  V.  (g)  The 
moment  tensor  Mo,  plotted  in  a  lower-hemisphere  projection.  The  location  of  the  piercing 
point  for  each  station  depends  on  the  station  azimuth,  epicentral  distance,  and  the  assumed 
layered  reference  model. 
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Figure  5.  Waveform  fits  for  four  nuclear  explosions,  (a)  BEXAR ,  (b)  HORNITOS,  (c) 
BRISTOL,  and  (d)  DIVIDER.  See  Figure  6  for  misfit  plots  on  the  lune. 
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Figure  6.  Lune  plots  for  four  nuclear  explosions. 

Misfit  plots  on  the  lune  for  four  example  nuclear  explosions  whose  best-fitting 
waveforms  are  shown  in  Figure  5.  The  color  scale  represents  the  variance  reduction,  from 
high  (blue)  at  the  best- fitting  moment  tensor,  to  low  (red);  the  min/max  values  are 
adjusted  in  each  subplot  to  show  the  variation  on  the  lune. 
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Station  BK.CMB.BH 
Mean  dist  380  km,  mean  az  287  deg 
Rayleigh:  mn  -12  md  -12  st  0.4 
Love:  mn  -02  md  0.5  st  22 


Station  US.MNV.BH 
Mean  dis*  216  km,  mean  az  308  deg 
Rayleigh:  mn  4.7  md  4.7  st  0.5 
Love:  mn  3.8  md  3.8  st  1.7 
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Figure  7.  Surface  wave  time  shifts  for  two  examples  stations  (CMB  and  MNV)  for  all  17 
nuclear  explosions  and  all  12  earthquakes  in  our  study. 


Figure  8.  The  influence  of  body  waves  on  lune  plots. 

The  influence  of  body  waves  and  surface  waves  on  moment  tensor  inversions  for  two 
example  events.  The  color  scale  is  adjusted  for  each  subplot  in  order  to  show  the 
variation  in  VR(A).  (a)  Little  Skull  Mountain:  body  waves  only;  waveform  fits  in  Figure 
S12.  (b)  Little  Skull  Mountain:  surface  waves  only,  (c)  Little  Skull  Mountain:  body 
waves  and  surface  waves,  (d)  HOYA:  body  waves  only,  (e)  HOYA:  surface  waves  only. 
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(f)  HOYA:  body  waves  and  surface  waves;  waveform  fits  in  Figure  3,  expanded  version 
of  uncertainty  summary  plot  in  Figure  4. 


(a) 


(«>) 


ISO 


ISO 


Full  moment  tensor  catalogs,  plotted  on  the  lune  (left)  and  on  the  vw  rectangle  (right) 
(Tape  and  Tape,  2012,  2015).  Beachballs  are  colored  by  magnitude,  (a)  32  earthquakes, 
explosions,  and  collapses  from  the  Nevada  Test  Site  (Ford  et  al.,  2009).  (b)  Same  events 
as  (a),  but  showing  our  solutions  in  this  study,  (c)  63  events  from  Uturuncu  volcano, 
Bolivia,  (d)  21  earthquakes  from  southern  Alaska. 
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